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ABSTRACT 

We introduce a multipolar scheme for describing the structure of stationary, axisym- 
metric, force-free black-hole magnetospheres in the "3+1" formalism. We focus here 
on Schwarzschild spacetime, giving a complete classification of the separable solutions 
of the stream equation. We show a transparent term-by-term analogy of our solutions 
with the familiar multipoles of flat-space electrodynamics. We discuss electrodynamic 
processes around disk- fed black holes in which our solutions find natural applications: 
(a) "interior" solutions in studies of the Blandford-Znajek process of extracting the 
hole's rotational energy, and of the formation of relativistic jets in active galactic 
nuclei and "microquasars" , and, (b) "exterior" solutions in studies of accretion disk 
dynamos, disk-driven winds and jets. On the strength of existing numerical studies, we 
argue that the poloidal field structures found here are also expected to hold with good 
accuracy for rotating black holes, except for maximum possible rotation rates. We 
show that the closed-loop exterior solutions found here are not in contradiction with 
the Macdonald-Thorne theorem, since these solutions, which diverge logarithmically 
on the hole's horizon H, apply only to those regions which exclude H. 

Key words: black hole physics - accretion disks - galaxies: active - galaxies: nuclei 
- galaxies: jets 



1 INTRODUCTION 

The study of black hole magnetospheres derives its impetus 
from our desire to understand the central powerhouse in ac- 
tive galactic nuclei (AGN) and double radio sources. In a 
seminal paper, Blandford and Znajek (1977, henceforth BZ) 
proposed that magnetic fields threading a Kerr black hole's 
horizon Ti can tap the rotational energy of the hole, trans- 
porting the energy in twin beams to power luminous dou- 
ble radio lobes. The structure of force-free magnetospheres 
surrounding Kerr black holes was discussed by BZ in a for- 
mal four-dimensional language, giving two specific examples 
of simple field structure: (a) split monopole and (b) non- 
separable paraboloid (see §2). 

In a pioneering effort to bridge the gap between the 
language of black-hole elctrodynamics and that of flat-space 
elctrodynamics with which astrophysicists are normally fa- 
miliar, Thorne and Macdonald (1982, henceforth TM) and 
Macdonald and Thorne (1982, henceforth MT) introduced 
the "3+1" formalism, wherein the spacetime was split up 
into three space directions and one, uniquely chosen, "uni- 
versal" time direction (see §2 for a brief sketch; a lucid expo- 
sition of the formalism is given in Thorne, Price & Macdon- 
ald 1986, henceforth TPM). Expressed in this way, equations 
of black-hole electrodynamics become remarkably similar to 
those of pulsar electrodynamics, and one can justifiably ex- 
pect to carry the accomplishment, expertise, and intuition 



in pulsar elctrodynamics over to black-hole electrodynamics. 

It is surprising, therefore, that a formalism for describ- 
ing the general field structure of black-hole magnetospheres 
has not been constructed so far, although the methods of 
the "3+1" formalism are ideally suited for this. Indeed, 
it seems that only four specific, simple field structures in 
Schwarzschild spacetime have ever been reported in the lit- 
erature, which appear to have been found ad hoc. Two of 
these are separable solutions of the stream equation (see 
§2): the monopole (and its trivial generalization, the "split 
monopole"; BZ; Macdonald 1984, henceforth M84) and the 
uniform field (Wald 1974; Hanni & Ruffini 1976; M84). 
The others are non-separable solutions: (a) asymptotically 
paraboloidal field lines (BZ), and, (b) asymptotically verti- 
cal field lines with specified slopes on the equatorial plane 
(Ghosh & Abramowicz 1997, henceforth GA). Apparently, 
the curved-spacetime analogue of our ail-too familiar and 
highy useful multipole description of electric and magnetic 
fields in flat-space electrodynamics has remained unnoticed 
so far. 

This paper introduces the first steps towards evolv- 
ing a general scheme for describing black-hole magneto- 
spheres in terms of their basic constituents, the multipo- 
lar "interior" and "exterior" solutions in the appropriate 
spacetime. We focus in this work on the field structure of 
stationary, axisymmetric, force-free magnetospheres around 
Schwarzschild black holes. We give a complete characteriza- 
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tion of separable solutions for the case in which the poloidal 
and toroidal fields are decoupled. We show that the angular 
parts of the stream functions which describe the magnetic 
field (see §2) are described by Gegenbauer polynomials. The 
radial parts of the stream functions divide themselves into 
two classes, interior and exterior solutions, in analogy with 
the nomenclature of flat-space electrodynamics. The former 
class is described by Jacobi polynomials, and the latter, by 
Jacobi functions of the second kind. We demonstrate a clear 
analogy between the basic constituents of the Schwarzschild 
solutions and those of flat-space electrodynamics, identify- 
ing Schwarzschild monopoles, dipoles, quadrupoles, and so 
on. We show that the existence of exterior solutions with 
closed field loops is not in conflict with a theorem proved by 
MT on the impossibility of closed field loops threading the 
horizon H, since these solutions, which diverge logarithmi- 
cally on Tt, are appplicable only to regions which exclude H, 
entirely in analogy with the behaviour of exterior solutions 
of flat-space electrodynamics with respect to the origin. In- 
deed, our inability to find any solutions which have closed 
loops and remain finite on TC is a confirmation of the MT 
theorem. 

Those studies of astrophysical processes involving black- 
hole electrodynamics which benefit particularly from this 
approach include (a) BZ-process studies with "realistic" 
field configurations suggested by recent numerical simula- 
tions (GA and references therein; particularly see Balbus & 
Hawley 1998; Brandenburg et al. 1995), (b) studies of mag- 
netic coupling of our exterior solutions to accretion disks in 
relation to dynamo action, jets and hydromagnetic winds 
(Konigl & Kartje 1994; Khanna & Camenzind 1992,1996; 
Blandford 1998), and, (c) studies of our interior solutions 
in relation to magnetohydrodynamic (MHD) models of jet 
formation and collimation in AGN and "microquasars" (Ca- 
menzind 1995; Fendt 1997, henceforth F97; Koide, Shibata 
& Kudo 1998, 1999, henceforth K98, K99; Eikenberry et 
al. 1998). Distinct field geometries relevant for these dif- 
ferent processes usually emerge naturally as low-order (see 
§5.4) examples of different classes of solutions identified in 
this work: we indicate these briefly, deferring detailed astro- 
physical applications to later works. Indeed, a main virtue 
of our multipolar approach to black-hole elctrodynamics is 
that a specific model field-configuration is immediately iden- 
tifiable as the natural choice for a given astrophysical prob- 
lem from symmetry considerations alone, without detailed 
computations. 

Calculations for rotating black holes will be given in 
a later work. We note here that, in a numerical study uti- 
lizing relaxation methods, M84 spun up Schwarzschild so- 
lutions corresponding to the first three of the four simple 
cases listed above, and showed that the effect of rotation 
on the poloidal field structure was very small in all cases, 
being barely pereceptible for < < 0.75, where is 
the hole's angular momentum per unit mass. If this result is 
generally true, it is of basic importance, since it would im- 
ply that, as far as the poloidal field structure is concerned, 
the hole's rotation is not just an inessential complication, 
but also an insignificant one except at the fastest possible 
rotation rates. For astrophysical applications, the poloidal 
structures classified in this paper thus appear to be ade- 
quate in all cases not involving black holes rotating near the 
limiting rate. Recent numerical simulations of jet formation 



by Koide and collaborators (K98; K99; Koide, Meier, Shi- 
bata & Kudo 1999, henceforth KMSK), using their general 
relativistic MHD code, support this conclusion: we return to 
this point in §5.4. In a recent numerical study of the stream 
equation utilizing finite-element techniques, F97 computed 
field structures (also see §5.4) in Kerr spacetime which join 
on smoothly to the asymptotic jet solutions given by Appl 
& Camenzind (1993). Finally, it must be clear on general 
grounds that the toroidal field strength is expected to de- 
pend significantly on the hole's rotation rate. We discuss 
these points again in §5.4, suggesting that the effects of the 
hole's rotation on magnetospheric structure may be of sec- 
ondary importance in studies of the BZ process and related 
issues, while they may be crucial in understanding collima- 
tion and stability of jets in AGN and microquasars. 

2 BLACK HOLE ELECTRODYNAMICS 

2.1 The "3+1" formulation 

Only the essentials of the "3+1" formulation of curved 
spacetime elctrodynamics are given here, since detailed ex- 
positions may be found in TM, MT, and TPM. The four- 
dimensional formulation of black hole electrodynamics can 
be reexpressed in terms of the familiar electric, magnetic, 
and current three-vectors and the charge density scalar of 
flat-space electrodynamics through the use of the following 
prescription. At each spacetime event, one chooses a fiducial 
reference frame by doing a "3+1" split-up of spacetime into 
three space directions and one, uniquely chosen, "univer- 
sal" time direction. One then decomposes the electromag- 
netic field tensor into the electric and magnetic field vectors 
in the usual manner, as also the four-current density into 
the current three-vector and the charge density, rewriting 
Maxwell's equations, the continuity equation for charge and 
current density, and the momentum equation in terms of 
these variables. 

In stationary black-hole electrodynamics, the natural 
choice for the fiducial frame is that of a so-called ZAMO 
(zero angular momentum observer; see Bardeen et al. 1973), 
who is at rest in the hole's stationary gravitational field, 
which, for rotating holes, implies further that the observer's 
angular velocity must be such that the observer's world 
lines appear orthogonal to a family of three-dimensional 
hypersurfaces of constant (Boyer-Lindquist) time (TPM). 
By mentally collapsing the entire collection of the above 
constant-time hypersurfaces into a single three-dimensional 
space, which we can regard as a sort of "absolute space", 
we can envisage electrodynamics (and other physics) as oc- 
curring in this absolute three-space, the universal time be- 
coming a parameter which demarks the evolution of mat- 
ter and fields in this absolute space. This is the absolute- 
space/universal-time picturization of the 3+1 formalism, 
whose analogy with Galilean relativity facilitates extensions 
of astrophysical intuition into black hole electrodynamics 
(TM). 

2.2 The stream equation 

In the above framework, the structure of a stationary, ax- 
isymmetric, force- free black hole magnetosphere is specified 
by a scalar function, if), of position, called the stream func- 
tion, and two functions of ip, namely the currrent potential, 
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J, and the angular velocity, Qf, of the magnetic field lines 
(M84). The poloidal and toroidal components of the mag- 
netic field are given by 



B r 



and 



Bt = 



Vtp x 4> 
2ircu 



(i) 



(2) 



respectively, with analogous expressions for electric fields, 
the toroidal electric field being zero in an axisymmetric sys- 
tem. Here, V is the gradient operator of the absolute three- 
space referred to above, i.e., the "space-space" part of the 
four-dimensional metric of the spacetime that we began with 
(M84). Further, K is the so-called lapse function or gravita- 
tional redshift factor, i.e., the rate of lapse of the ZAMO's 
proper time relative to that of the universal time referred 
to above: it is the "time-time" part of the four-dimensional 
metric. Finally, vj = rsin# is the cylindrical radius, and <f> 
is a unit vector in the azimuthal direction. 

The name "stream function" was (apparently) coined 
by Newtonian analogy, since the poloidal magnetic field 
points everywhere along streamlines of tp (i.e., poloidal lines 
of constant tp): thus, poloidal field lines are contours of con- 
stant tp [see eq.(l)]. Uses of this description for magnetic 
fields are well-known in flat-space electrodynamics : see, 
e.g., Michel's (1973, henceforth M73) work on pulsar elec- 
trodynamics, in which tp was called a field-line label. Fur- 
ther analogies with flat-space elctrodynamics may be made: 
tp is proportional to the toroidal component of the vector 
potential (which, in turn, is related to the magnetic flux), 
and I(n),9o) is the total current through the azimuthal loop 
(r = r , 6 = 6q). The stream function satisfies a second- 
order elliptic differential equation called the stream equation 
(MT), and also, sometimes, the Grad-Shafranov or trans- 
field equation (Beskin 1997, henceforth B97). The equation 



K 2 



n_F — Ul d^F /„ , n? 

H -(ViA) 



16n 2 r dl 



Kro 2 dtp 



I-rr = 0. 



(3) 



Here, lo is the angular velocity of the ZAMO rest frame 
relative to absolute space, i.e., the "time-space" part of the 
four-dimensional metric referred to above. 

For Schwarzschild black holes, £If = = u>, so that the 
stream equation simplifies considerably. The lapse function 
is given in this case by 



N = Wl- 



2M 



(4) 



where M is the hole's mass. (We shall use the natural units, 
G = 1 = c throughout this paper.) In this work, we study 
those black-hole magnetospheres in which the current poten- 
tial I(tp) is of such form that the last term on the left-hand 
side of eq.(3) vanishes, so that the poloidal and toroidal 
fields are, in some sense, decoupled. (For a different choice 
of I(tp), useful for modeling cylindrically collimated jets, see 
Appl & Camenzind 1993 and F97). Two types of (astrophys- 
ically) important situations in which this condition holds 



have been described in the literature. The first type occurs 
when the black hole is immersed in a vacuum: all analytic 
Schwarzschild solutions reported by Wald 1974, Hanni & 
Ruffini 1976, BZ, and M84 fall in this category. The second 
type corresponds to the situation described by GA, in which 
/ = constant = Io, so that Bt oc (see eq. [2]). In this 

case, the toroidal field scales asymptotically as B<j, ~ zu^ 1 
far from the horizon H and diverges, of course, as (MT) 
as H is approached. This asymptotic behavior is reminis- 
cent of some of the Newtonian model solutions of Blandford 
(1976, henceforth B76), and is relevant for electrodynamics 
of accretion disk-fed black holes, as emphasized by GA. In 
fact, these two types of situations are, to our knowledge, 
the only ones in which analytic Schwarzschild solutions of 
field structures around black holes have been reported so 
far. In these circumstances, the stream equation reduces to 
the linear differential equation 



{>}- 



0. 



(5) 



2.3 Separable solutions 



With the aid of equation (4) and the explicit form of V for 
the Schwarzschild metric, equation (5) becomes 



dr 



(1 



2M dtp 
~~r>~dr 



86 



1 dtp 



0. 



(6) 



In this work, we are principally concerned with separable so- 
lutions of equation (6), which offer a rich variety of multipo- 
lar structure. Only two non-separable solutions of equation 
(6) have been reported in the literature (BZ, M84, GA): we 
refer to these at appropriate places. 

Separable solutions of equation (6) are of the form 



tp(r,6) = R(r)@{e) , 



(7) 



where the radial part, R, and the angular part, 6, satisfy 
equations which are obtained by combining equations (6) 
and (7): 



d6 
and 
d_ 

dr 



1 d@ 



sin (9 dO 



_ 2M_ dR 
r dr . 



-A- 



e 



(8) 



(9) 



Here, A is the separation constant. 



3 SOLUTIONS OF THE STREAM EQUATION 
3.1 Lowest-order solutions 

Solutions of equations (8) and (9) in the special case A — 
constitute the lowest-order separable solutions of the stream 
equation. In the next subsection, we quantify the idea of 
order in terms of the value of A. 

The lowest-order solutions are obtained immediately on 
setting A — in eqs. (8) and (9). These are 



O = a cos 6 + b , 
and 

R = c[r + 2M ln(r - 2M )] + d . 



(10) 



(11) 
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Here, a, b, c, d are constants. The nature of the solutions is 
discussed in detail in §4. Note that the special case c = 0, 
the Schwarzschild monopole, was one of the first solutions 
to be found. Note further that the special case a + b = 
0, d — is that of the separable Schwarzschild paraboloid, 
the non-separable analogue of which has been reported (BZ; 
M84) and is discussed later. Finally, note the logarithmic 
singularity at r = 2M in eq.(ll), which is discussed later. 



3.2 Solutions of general order 

We now give a complete characterization of general solutions 
of equations (8) and (9), to be labelled by the ordinal num- 
ber, m, which is defined in terms of the separation constant 
as 



A = m(m + 1) . 



(12) 



With the aid of equation (12), the equation for the angular 
function, 0, becomes 



(1-a; 2 )— + m(m+ 1)6 = 0. 



(13) 



Here, x = cosd. Similarly, the equation for the radial func- 
tion, R, becomes 

(l-z 2 )^4~^ + rn(rn + l)R = 0. (14) 
dz 2 dz 

Here, z=s.-l. 



3.3 The angular function 

Solutions of equation (13) are closely related to the Gegen- 
bauer polynomials (see, e.g., Szego 1939, henceforth SZ, 
Abramowitz & Stegun 1972, henceforth AS; Gradshteyn & 
Ryzhik 1990, henceforth GR), the angular function of order 
m being given by 



e m (x) = (i - x 2 )c l *Ux) 



(-) 

where C m 2 _ x are Gegenbauer plynomials. Explicit forms for 
(-) 

C m 2 _i can be obtained with the aid of recursion relations, 

(-) (-) 

starting with expressions for C 2 andC x 2 (SZ). We list be- 
low the first few. Note that we shall henceforth suppress the 

(3) 

superscript in CJ 2 _ 1 , since its value is the same throughout 
this work. 

m = 1 Co (a') = 1 , 



(15) 



3.4 The radial function 

Solutions of equation (14), which are closely related to the 
Jacobi polynomials and Jacobi functions of the second kind 
(see SZ for an excellent review), are given by (AS, GR, SZ) 

m-l( Z ) 



where "P^l^are Jacobi polynomials and 2^{are Jacobi 
functions of the second kind. The radial function of order 
m is therefore given by 

' 'n(°' 2 )(S. 1) 

' m-lV M > 



l(0.2) . 



R m (r) 



n(°' 2 ) is. i\ 



(18) 



The Jacobi polynomials and the Jacobi functions of the sec- 
ond kind represent two distinct classes of radial solutions, 
whose physical natures are evident from their asymptotic 
behaviors at large r. For r S> M , their scalings are (AS, SZ) 



V 



(0,2) 



O (0 - 2) ~ r - m -' 2 



so that the radial function scales as 



R, 



.W = { 



(19) 



(20) 



in the two cases, which is immediately reminiscent of the 
behaviors of the "interior" and "exterior" solutions of flat- 
space elctrodynamics. Indeed, since r > M is the flat-space 
limit of Schwarzschild results, the analogy is quite exact. 
Therefore, we henceforth refer to the solutions bearing the 
Jacobi polynomials as the interior solutions and those bear- 
ing the Jacobi functions of the second kind as the exterior 
solutions. Since the two types of solutions have different 
physical significances and mathematical properties, we treat 
them separately. 

3.4-1 Interior solutions 

Explicit forms for Jacobi polynomials can be obtained from 
the above references: we list below the first few. Note that 
henceforth we suppress the two superscripts in V^fl, as 
their values are the same throughout this work. 

m = l Vo(z) = l, 

m = 2 P 1 (z) = 2z-1, 

1 C C 1 

m = 3 V^z) = ^-\z-\. (21) 

For our purposes, it is useful to scale the radial function 
Rm(r) in terms of its asymptotic value given in equation 
(20). We combine eqs.(17) and (21), and obtain 



R m (r) = r m+1 g m ( — ) 



(22) 



■2 Ci(z)=3x, 



m = 3 Ci(£) 



15 2 
~2 X 



(16) 



Detailed angular structure of the solutions is discussed in 
§4. Note that the m = 1 solution has dipolar symmetry, 

01 ~ sin 2 6, the m = 2 solution has quadrupolar symmetry, 

02 ~ sin 2 8 cos 6, and so on. This is in keeping with our 
experience in flat-space electrodynamics. 



where Q m {^r) is a polynomial which describes the general- 
relativistic effects. The first few of these polynomials are 

rri — 1 Si{x) = 1 , 



m — 2 Qi{x) = 1 - -x , 

8 8 

m = 3 Qz{x) = 1 - -x + - x . 

o O 



(23) 



The structure of the interior solutions is discussed in §4. 
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3.4-2 Exterior solutions 

As the explicit forms for Jacobi functions of the second kind 
are not easily found in the standard references, we describe 
them here at the level of detail adequate for our purposes. 
We introduce a prescription for obtaining a general func- 
tion of this type, Q^ ,l3 \ which, to our knowledge, is not 
found in the standard mathematical references. We suggest 
that explicit forms for Q^'^'can be found most conveniently 
by the following prescription, through the use of recursion 
relations, starting with a suitable (integral or differential) 
representation for Q^'^. 

For arbitrary values of a, /3 satisfying the constraint 
a + (3 + 1 > 0, a useful integral representation (SZ) for 
is 



Q^\z) = -2' 



t+/3 i> + i)r(/3 + i) 
r{a + p + i) 



f 



dt 



(24) 



For positive integral values of a and /3, as in our case, a 
differential representation, which leads to much faster cal- 
culation of explicit expressions for Qq*'^, can be obtained 
from equation (24) . The representation is 



d a of f 1 t z-a ] 
da a db' 3 \ a + b U z + b } 



(-) 



/3 + l 9 a+/3 



r(a + /3 + l) 



(25) 



and we have not found it in any standard reference. To prove 
equation (25), consider the identity 



f 

J OO 



dt 



In - 



(t-a)(t + b) a + b z + b 



(26) 



Differentiation of both sides of equation (26) a times with 
respect to a and (3 times with respect to b, and comparison 
with equation (24), leads to equation (25). 

Recursion relations from which o!*n'^ with n > can 
be obtained are as follows (SZ). The relation used for calcu- 
lating Q[ a '^ is a special one: 



Q { t P) (*) = + P + 2)z + a - /3\Q£" t "{z) - 2 



(«,/3), 



■>a+/3-l 



(27) 



a + (3 + 2 r(a + l)r(/3 + l) 
{z-l)<*{z+\y T(a + /3 + 2) 

For n > 2, Q^ p) are obtained by repeated application of 
the general relation which holds for n > 2, namely 



2n(n + a + /3)(2n + a + f3 - 2)Q { n °" p> ( 



(2n+a+f3-l)[(2n+a+l3)(2n+a+/3-2)z+a'-0']Q ( ^ ,3 1 \z) 



-2(n + a - l)(n + (3- l)(2n + a + /^Q^^) 



(28) 



We now set a = and /3 = 2 for our problem, and cal- 
culate Jacobi functions of interest from equations (25), (27) 



and (28). We give the first few below, henceforth suppress- 
ing the two superscripts in Q^l\ , since their values remain 
the same throughout this work: 

z + 2 

2 - T ~ ( 2 +l) 2 ' 



i r\ ( \ 1 i z + 1 
m = 1 Qo(z) = - In 



m = 2 Qi(x) = («- 



In 



2+1 



Z - 1 



2« 2 + 3z + § 

(2 + 1) 2 



m = 3 62(2) = (-^-z 2 



1.. 2 + 1 

8 )ln I^T 



if2 3 +52 2 



(29) 



The structure of the exterior solutions is discussed 
in §4. Note the appearance of a logarithmic singularity, 
~ ln(r — 2M) (z in eqs.[29] is given by z — -p- — 1; see 
above), in these solutions at the horizon fi; the singularity 
also appeared in one of the lowest-order solutions discussed 
above. This singularity is reminiscent of the one which oc- 
curs in the well-known flat-space exterior solutions (~ r~ m ) 
at r = 0, but it is weaker. Of course, such a singularity does 
not imply the invalidity of exterior solutions for black hole 
electrodynamics, any more than the above power-law singu- 
larity implies the invalidity of exterior solutions in flat-space 
electrodynamics: it simply means that these exterior solu- 
tions are valid in regions of space which exclude TC, just as 
the regions of validity of flat-space exterior solutions nec- 
essarily exclude the origin. We come back to this point in 
§4 and §5. Finally, we note that the power-law singularity 
displayed by the Jacobi functions of equation (29) at r = 
(i.e., z + 1 = 0) has no physical significance, since these 
solutions do not apply for r < 2M. 



4 STRUCTURE OF THE SOLUTIONS 
4.1 Lowest-order solutions 

The solution corresponding to c = in equation (11), i.e., 



ip ~ a cos 9 + b , 



(30) 



with a and b constant, is the Schwarzschild monopole, the 
field lines of which are shown in Figure 1. 

This solution has been discussed previously by BZ and 
M84, and its flat-space version (identical in form to eq.[30]) 
by M73. Imposition of the solenoidal condition on B, which 
must be satisfied in all astrophysical applications, is achieved 
by the device of assigning opposite polarities to the radial 
magnetic field in the two hemispheres (BZ, M84), i.e., a 
so-called "split monopole" configuration, maintained by a 
toroidal current in the accretion disk lying in the equatorial 
plane. Results of spinning up the split monopole field have 
been described in various approximations by BZ and M84, 
and their flat-space analogue is an exact solution given by 
M73. 

The solution corresponding to a + b = 0, d — 0, i.e., 



1> 



+ 21n(r-2M)][l-cos0], 



(31) 



represents the separable Schwarzschild paraboloid (our 
nomenclature comes from the observation that, in eq.[31], 
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Figure 1. Poloidal field lines of the Schwarzschild monopole. The 
horizon Tt is shown by the thick line. All lengths are in natural 
units, as indicated. 



Figure 2. Same as Figure 1, but for the separable Schwarzschild 
paraboloid. 



the sole difference from a familiar flat-space paraboloid 
is the presence of a logarithmic term characteristic of 
Schwarzschild spacetime) , the field lines of which are shown 
in Figure 2. 

At sufficiently large radii, this configuration reduces 
to the paraboloidal field-lines of flat-space electrodynamics, 
which have been studied at length by B76. At small radii, 
the separable Schwarzschild paraboloid has one curious rel- 
ativistic feature: close to Ti, the field lines change direction, 
so that tp, which is a measure of the magnetic flux (see §2), 
changes sign. For the specific stream function given in equa- 
tion (3f), ij) changes sign at r m 2.314M. Note that this 
asymptotically-paraboloidal separable solution has similari- 
ties to the corresponding non-separable solution (BZ, M84), 
which we discuss in §5, but there are also clear differences. 

4.2 General solutions 

We now combine the angular and radial parts of general 
solutions for the stream function described in §3, and discuss 
the structure of the resultant field configurations. 



4- 



Interior solutions 



The stream function for m = 1, namely, 
rpi ~ r 2 sin 2 9 , 



(32) 



represents a uniform magnetic field pointing in the z- 
direction, as shown in Figure 3. 

This appears to be the most thoroughly studied (Wald 
f974; Hanni & Ruffini f976; M84) field configuration in 
Schwarzschild spacetime, because of its simplicity. The con- 
figuration is identical to that obtained in flat-space electro- 
dynamics corresponding to the m = 1 term of a multipolar 
decomposition of the scalar potential, i.e., $1 = rPi(cos(9) 
(here Pi is a Legendre polynomial). 

The stream function for m = 2, namely, 

V>2 ~ (r 3 - ^Mr 2 ) sin 2 9 cos 9, (33) 

represents the interior solution of the next higher order, the 
field lines of which are shown in Fig.4. 

This configuration is very similar, but not identical, to 
that obtained in flat-space electrodynamics corresponding 
to the m — 2 term of a multipolar decomposition of the 
scalar potential, i.e., $2 = r 2 P2(cos#). 

The structure of higher-order solutions may be read- 
ily obtained in a similar manner. It is this term-by-term 
analogy with well-known configurations in flat-space elec- 
trodynamics that makes our scheme of classifying the field- 
configurations in curved spacetime particularly transparent, 
and underscores the power of the 3+1 formalism. 

4-2.2 Exterior solutions 
The stream function for m = 1, namely, 
r 



if) 



2 r 



2M 



-rM-Af 



Figure 3. Same as Figure 1, but for the uniform field in 
Schwarzschild spacetime. 



Figure 4. Same as Figure 1, but for the m = 2 interior solution. 

represents a solution which we name the Schwarzschild 
dipolc. The field lines, shown in Figure 5, readily explain 
the terminology, demonstrating that this solution is the 
Schwarzschild analogue of the familiar flat-space dipole field. 

The analogy is quite exact since, on taking the limit r ^> 
M in eq.(34), we find that ipi ~ i sin 2 9, which describes the 
field lines of a flat-space dipole. 

Similarly, the stream function for m = 2, namely, 



V>2 ~ [(r 3 - I Mr 2 ) In — 



M 3 



2M 



2r 2 M + rM 2 + 



(35) 



sin 2 9 , 



(34) 



represents the Schwarzschild quadrupole, the field lines of 
which are displayed in Figure 6. 

The quadrupolar symmetry is evident. The reader can 
verify the exactness of the flat-space analogy by extracting 
the r > M limit of equation (35), and noting that, in this 
limit, ip2 = const is precisely the equation for the field lines 
of a flat-space quadrupole. 

The higher exterior multipoles (octupole and so on) may 
be readily obtained in a similar manner. For the exterior 
solutions, the exactness of the analogy with the multipoles 
of flat-space electrodynamics makes the classification scheme 
completetly natural and transparent, and re-emphasizes the 
virtue of the 3+1 formalism. 



5 DISCUSSION 

5.1 Astrophysical context 

The classes of solutions detailed in §§3-4 provide natural 
generic choices for appropriate classes of black-hole electro- 
dynamics problems involving accretion disks. For studies of 
field structure and transport of energy or angular momen- 
tum close to the hole's horizon, our lowest-order solutions 
and interior solutions (BZ; M84; B97; KMSK) are clearly the 
best choice (also see §5.4). The uniform- field configuration, 
corresponding to m — 1, describes the essential astrophysics 
very close to H, as has been intuitively realized long ago 
(see TPM, particularly their Fig. 36). This is a result of the 
physical mechanisms described in MT and TPM by which 
the black hole "cleans" the magnetic field threading its hori- 
zon of details, kinks, and closed loops (also see §5.3). At 
slightly larger distances from H (~ TZh, say, where TZn is 
the characteristic scale size of TL), further details of field 
structure begin to appear even under the most ideal circum- 
stances, due to the contribution of the currents circulating 
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Figure 5. Same as Figure 1, but for the Schwarzschild dipolc. 



Figure 6. Same as Figure 1, but for Schwarzschild quadrupolc. 

at the inner edge of the accretion disk. These are described 
by interior solutions of higher order, m = 2,3,...., whose 
relative contributions determine the shape of the beam or 
jet of relativistic particles transporting energy to the radio 
lobes. Similar arguments apparently also hold for MHD jets 
in AGN and quasars, although the stream equation there is 
somewhat different (B97 and references therein), and so is 
the structure of the interior solutions. 

On the other hand, for describing the field structure on 
and around the accretion disk at larger distances (S> Tin) 
from H, our exterior solutions are the natural choice. They 
describe, term by term, the multipole moments of the cur- 
rents circulating in the disk, seen as an external observer 
comes closer and closer to the hole. The m = 1 dipole is 
thus a first estimate of the total effect of the disk currents, 
as seen by a distant observer. This decomposition proves 
particularly useful for analyzing the relative effects of cur- 
rents circulating in the disk and on the hole's horizon. Astro- 
physical applications of these field structures to disk-driven 
jets and MHD winds, and to dynamo-action in disks are 
well-known (Konigl & Kartje 1994; Khanna & Camenzind 
1992,1996; Camenzind 1995). 

Construction of model global field-structures relevant 
to almost all black-hole-accretion-disk problems involves the 
matching of interior and exterior solutions at a suitable in- 
termediate radius according to a prescription appropriate 
for the particular problem. This has direct relevance to the 
"impedance matching" considerations long known to be of 
crucial for studies of BZ processes and related particle ac- 
celeration issues (Lovelace, MacAuslan & Burns 1979; MT; 
TPM). In fact, techniques borrowed from matching proce- 
dures in flat-space electrodynamics may provide the most 
rigorous formulation of the idea of impedance matching be- 
tween the black hole source and the distant astrophysical 
load: we shall describe this elsewhere. 

5.2 Non-separable solutions 

Separable solutions do not exhaust all possibilities for solu- 
tions of equation (5) . One particular non-separable solution 
of astrophysical interest is that found by BZ and discussed 
by M84. It is given by 

ip ~ [(r - 2M)(1 - cos (9) - 2M(1 + cos (9) ln(l + cos (9)] ,(36) 

and its field lines are asymptotically paraboloidal. This so- 
lution was used in the original BZ study of the extraction 
of rotational energy from black holes. Equation (37) is one 
of the lowest-order, non-separable solutions of equation (5) , 
analogous to the lowest-order separable paraboloid solution 
described by equation (31). At large radii, both of these 
solutions reduce to the familiar, exact paraboloids of flat- 
space electrodynamics, ip nr ~ r(l — cosO), astrophysical ap- 
plications of which to accretion- disk electrodynamics and 
jet-formation have been discussed at length by B76. 



However, there is nothing fundamental about the 
paraboloidal shape in the astrophysical context, except that 
field lines do need to be asymptotically vertical if electro- 
magnetic energy is to be beamed along the black hole's ro- 
tation axis, as in the BZ process. In fact, the asymptoti- 
cally paraboloidal shapes discussed above have the disad- 
vantage that they have an inflexible shape on the equato- 
rial plane (the ratio B m jB z is always unity on this plane 
in the Newtonian limit), so that it is sometimes impossi- 
ble to satisfy boundary conditions of astrophysical interest 
on the accretion-disk surface. A good example is discussed 
in GA: recent numerical simulations of dynamo-generated 
magnetic fields in accretion disks suggest B m /B z ^ 2 to 3 
(GA and references therein), which paraboloidal field lines 
cannot mimic. To achieve this, GA found a solution 

i> ~ [(r - 2M)(1 - cos 6) - r cos 6 

-2M(l + cos0)ln(l + cos0)], (37) 

which was adequate for this kind of modelling. It is clear that 
the solution given by equation (37) is a linear superposition 
of those given by eqs. (36) and (30), and therefore must nec- 
essarily be a solution of the linear equation (5). Equations 
(36) and (37) appear to be the only astrophysically inter- 
esting non-separable solutions of equation (5) known at this 
time. 

There is no difficulty in obtaining further non-separable 
solutions by superposing equation (36) with other classes 
of separable solutions found in this work, but this is not of 
basic importance. A systematic method of classifying higher- 
order non-separable solutions would be more important both 
in principle and in the astrophysical context: this will be 
discussed elsewhere. 

5.3 The MT theorem 

MT proved an intersting theorem (henceforth referred to as 
the MT theorem) which states that closed loops of magnetic 
field threading the horizon TL cannot exist in stationary, ax- 
isymmetric, force-free magnetospheres. Therefore, a natural 
question which may arise at this point is: how do the ex- 
ternal solutions found here, all of which have closed loops 
extending to TL, stand vis a vis the MT theorem? 

To answer this question, consider first the proof of 
the MT theorem. Since we are concerned here only with 
Schwarzschild black holes, we shall be content with the proof 
for this case, which is particularly simple. Since equation (5) 
implies that the vector ^Vtp is solenoidal, we can convert 
it , by Gauss's theorem, into the condition for vanishing of 
its surface integral over a closed surface C which consists of 
a closed loop and the region between its footpoints on TL: 

/ — V^-dS = 0. (38) 
Jc m 

Here, dS is an element of surface area, its direction be- 
ing along the outward normal. Now, MT considered stream 
functions which (and whose derivatives) remained finite on 
TL, so that their next argument was that, since the lapse 
function ft vanishes on TL by definition (see, e.g., Bardeen 
1973), so does that part of the integral which comes from 
TL, implying that the integral over the closed field loop itself 
vanishes. But, since the field lines are contours of constant tp 
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(see §2.2), Vy> is everywhere parallel to dS, so that the inte- 
grand has the same sign all over the closed loop. It follows, 
therefore, that the integral can vanish only if the integrand 
vanishes everywhere, which means that there can be no field 
or field loop. 

Consider now the results of applying the MT theorem 
to our closed-loop external solutions. Since the stream func- 
tions ip for these solutions diverge like In H as TL is ap- 
proached, as explained earlier (see §3.4.2 and eq.[4]), Vtp 
diverges as Thus the integrand, -t-VV>, is finite over 

TL, and so is its contribution to the integral in eq. (38). De- 
pending on the details of the field geometry in the close 
neighborhood of TL, it may, therefore, be possible to satisfy 
eq.(38) without requiring that the field be zero everywhere 
on the closed loop. Thus, the MT theorem is not in contra- 
diction with closed-loop solutions that diverge on TL. The 
situation is closely analogous to that for exterior solutions 
in flat-space electrodynamics: since these solutions, which 
have closed field loops, diverge at the origin (r = 0), inte- 
grals similar to that in eq.(38) (or related integrals involving 
the scalar potential) have to be evaluated carefully in the 
neighbourhood of the origin, and can have finite contribu- 
tions from there. 

The above discussion does not imply, of course, that 
the exterior solutions found here can extend upto TL. This is 
ruled out simply because they diverge on TL and are, there- 
fore, physically inadmissible there. As we have indicated ear- 
lier, their domain of validity is a region of space which ex- 
cludes TL, as the name exterior solution implies. Examples 
of this abound in flat-space electrodynamics, and it is not 
difficult to envisage model problems in black-hole electrody- 
namics where an inner region containing TL admits of inte- 
rior solutions, surrounded by an outer region which admits 
of exterior solutions: we have already indicated above that 
matching these solutions may be of considerable astrophys- 
ical importance. 

Finally, the MT theorem only excludes the existence of 
closed-loop field configurations which remain finite on TL. We 
have shown above that all separable Schwarzschild solutions 
with closed loops diverge on TL. Since the two non-separable 
solutions known so far (see above) do not have closed loops, 
it is correct to say that all known Schwarzschild solutions 
with closed loops diverge on TL. This is thus a confirma- 
tion of the MT theorem, showing that the configurations 
that MT showed to be untenable are, indeed, not found am- 
nong the known solutions. Again, we find here a close anal- 
ogy with multipolar expansion in flat-space electrodynam- 
ics, where all of the closed-loop multipole solutions (dipole, 
quadrupole, etc) diverge at the origin. 

5.4 Concluding remarks 

Our confidence in the relevance of the poloidal Schwarzschild 
solutions, found in this work, to electrodynamics around ro- 
tating black holes stems largely from the insignificant in- 
fluence of rotation, except at the highest possible rates, on 
poloidal field structures found in all cases investigated so far 
(M84; K98; K99; KMSK; also see below). A remaining con- 
cern is about the role of the toroidal field, whose strength is 
expected to depend significantly on the hole's rotation rate 
(MT). While poloidal and toroidal fields are effectively de- 
coupled, as we argued in §2, in most problems of interest 



in black-hole electrodynamics, so that the higher toroidal 
fields generated by increasing rotation rates do not react 
back on the poloidal fields, their influences on the essential 
astrophysical phenomena, e.g., the BZ process, need to be 
clarified. The key point here is the effect of the well-known 
boundary conditions on TL (MT and references therein). As 
we approach TL, the poloidal field becomes entirely normal, 
B±, to TL, while the toroidal field Bt diverges as so 
that Bh = NBt remains finite. The strength of the BZ 
process scales with the magnetic stress B±Bh, but the 'de- 
generacy' condition (E.B = 0) outside TL fixes the ratio 
B±/Bh in terms of the hole's rotation rate (GA and refer- 
ences therein), so that the BZ power output finally scales 
with B\ alone. Thus, any direct influence of the magnitude 
of the toroidal field on the strength of this particular as- 
trophysical process ultimately drops out, and a sufficiently 
accurate estimate of the poloidal field close to TL is all one 
needs to calculate the strength of the BZ process and re- 
lated quantities for a black hole with a specified rotation 
rate. This is why our poloidal field calculations are expected 
to be of considerable practical value in BZ-process calcula- 
tions relevant to double radio sources, AGN, and quasars, 
as also to the poloidal field structures associated with rela- 
tivists jets in AGN and microquasars. However, there are 
other astrophysical phenomena that do depend crucially on 
the strength of Bt away from the horizon, such as the col- 
limation and stability of the above jets (Appl & Camenzind 
1992; K99; KMSK). We expect, therefore, that adequate in- 
vestigations of these phenomena will require reliable calcu- 
lations of the toroidal field structure of the magnetospheres 
of rotating black holes, which is beyond the scope of this 
work. 

After the original version of this paper was submit- 
ted, we became aware of an excellent contaporaneous re- 
view (B97) on axisymmetric stationary flows in astrophys- 
ical objects, in which a model black-hole magnetosphere 
is constructed by expanding the nonrelativistic M73 solu- 
tion referred to earlier, and introducing relativistic effects 
into each term separately. The final result, given upto two 
terms, seems to be similar to a combination of our interior 
solutions with m — 1,3; the m = 2 term vanishes for this 
particular solution considered by B97 because of its sym- 
metry properties. We also note the similarity between the 
poloidal structure of our separable Schwarzschild paraboloid 
(Fig. 2), and that of the numerical computation in F97 of 
collimated jet magnetospheres in Kerr spacetime (Fig. 3 of 
F97): this may well be indicative of the secondary impor- 
tance of the hole's rotation on magnetospheric structure 
that we have suggested in this paper. Finally, the results 
of recent simulations of jet formation by Koide and col- 
laborators (K98; K99; KMSK) with the aid of their gen- 
eral relativistic MHD code clearly support our conclusions. 
The poloidal field structures around rapidly-rotating Kerr 
holes (KMSK) are found by these authors to be very simi- 
lar to those around Schwarzschild holes (K98), particularly 
for co-rotating Keplerian disks, confirming a major argu- 
ment used in this paper. Further, these authors start their 
simulations with a uniform poloidal magnetic field (our in- 
terior solution with m = 1, as shown in our Fig. 3), and 
obtain poloidal structures at the end of their runs which 
appear remarkably similar to our separable Schwarzschild 
paraboloid (a lowest-order solution shown in our Fig. 2), for 
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both Schwarzschild holes (Fig. 2 of K98) and rapidly-rotating 
Kerr holes (Fig. 1(d) of KMSK) accreting from co-rotating 
Keplerian disks. This gives us much confidence in the rele- 
vance of the basic "building block" poloidal field structures 
found in this work to essential astrophysical processes near 
black holes: future simulations should fully exploit our mul- 
tipolar expansion scheme. 
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